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ABSTRACT 


The solar wind is observed to display high speeds in high heliolatitude coronal holes and low speeds near the 
ecliptic plane. The heliospheric magnetic field associated with the solar wind plays a very important role in the 
transport and modulation of charged energetic particles, including galactic cosmic rays (GCRs) and solar energetic 
particles (SEPs), in the three-dimensional heliosphere. In previous studies, a constant solar wind speed, which 
is independent of heliolatitude, was assumed and commonly used in modulation modeling of cosmic rays. In 
this work, we investigate the realistic latitudinally dependent solar wind speed and utilize the theoretical models 
in hyperbolic and piecewise polynomial forms to explore the important effects on the heliospheric magnetic 
field and the diffusion coefficients (parallel, perpendicular, and drift scale) of cosmic rays in the presence of 
adiabatic focusing. Comparisons of the diffusion coefficients derived from standard Parker field and modified 
magnetic fields are presented. Since the structures and properties of different solar wind sources (coronal streamer 
belt, polar coronal hole, and transition region between them) differ from each other in essence, we suggest that 
the latitudinally dependent solar wind speed and the corresponding heliospheric magnetic field and diffusion 
coefficients with adiabatic focusing should be employed in the global modeling studies of GCRs and SEPs in 
the heliosphere. 
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1. INTRODUCTION 


The discovery of solar wind opened the door for space physics 
and symbolized the start of the fruitful space exploration age. 
A large number of authors have contributed to this important 
subject since then. As a seminal work, Parker (1958) first 
predicted radial, supersonic solar wind with typical speeds 
of hundreds of kms~!, which was soon proven by Mariner 
2 and other spacecraft. Through several decades of efforts in 
the community of space science with spacecraft observations 
and theoretical investigations, we have achieved a much better 
understanding of the solar wind, either spatially in the global 
heliosphere or temporally with the whole solar cycle. It has been 
widely agreed that fast solar wind with a characteristic velocity 
of 800km s7! originates from polar coronal holes (e.g., Zirker 
1977; McComas et al. 1998a), and generally, slow solar wind 
with a typical speed of 400km s7! is associated with coronal 
streamer belts near the ecliptic plane (e.g., Hundhausen 1977; 
Gosling et al. 1981; McComas et al. 1998b). 

As we know, the heliospheric magnetic field is “frozen into” 
the solar wind and travels along with it to form the Archimedes 
spirals (Parker 1958). Some properties of the Parker magnetic 
field model have been supported by observations (e.g., Smith 
et al. 1986; Forsyth et al. 1996), while some deviations of the 
magnetic field are also known (e.g., Schwadron & McComas 
2005). In the heliosphere, the average magnetic field is not 
constant with the heliocentric distance, and this spatially varying 
average magnetic field gives rise to the adiabatic focusing effect 
of charged energetic particles (e.g., Roelof 1969; Earl 1976; 
Schlickeiser & Shalchi 2008; Schlickeiser & Jenko 2010; He 
& Wan 2012a, 2012b; He & Schlickeiser 2014). The focusing 


length is defined by 


L'@) = he 0B(z) 
B(z) 0z 

where B(z) is the guide magnetic field with direction z. 
The adiabatic focusing effect should be considered when we 
investigate the diffusion coefficients and transport processes of 
charged particles, such as galactic cosmic rays (GCRs) and 
solar energetic particles (SEPs), particularly in the scenario of 
the inner heliosphere. 

The diffusion behavior of charged energetic particles in a 
large-scale turbulent magnetic field remains one of the funda- 
mental problems of long-standing importance in space physics, 
astrophysics, and plasma physics. Basically, achieving a better 
understanding of this problem is essential for studying the mod- 
ulation of GCRs, the propagation of SEPs, the diffusive shock 
acceleration, and laboratory plasmas and controlled fusion. In 
general, the diffusion tensor of charged particles in a turbulent 
magnetic field can be separated into a symmetric part consist- 
ing of parallel diffusion and perpendicular diffusion along and 
across the mean magnetic field, respectively, and an asymmet- 
ric part, i.e., the drift effect. The parallel diffusion plays an 
evidently important role in the transport of cosmic rays, so it 
has been extensively studied (e.g., Jokipii 1966; Bieber et al. 
1994; Schlickeiser & Shalchi 2008; He & Qin 2011; He & Wan 
2012a). The effect of perpendicular diffusion has also become 
increasingly well known in the research fields of cosmic rays 
and SEPs, and it attracts a number of authors (e.g., Dwyer et al. 
1997; Giacalone & Jokipii 1999; He et al. 2011; He & Wan 
2012b). Recently, He et al. (2011) found that the perpendicu- 
lar diffusion plays a very important role in the propagation of 


(1) 
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SEPs, particularly when an observer is not directly connected 
to the acceleration regions by the interplanetary magnetic field 
lines. In light of the direct method (He & Wan 2012a, 2012b), 
He & Wan (2013) theoretically investigated the rigidity depen- 
dence of the parallel and perpendicular mean free paths, in an 
event-to-event view, and basically reproduced the observational 
results of several typical SEP events including a “counterex- 
ample” presented in previous well-known works (e.g., Palmer 
1982; Bieber et al. 1994). They also reported that when tak- 
ing a series of SEP events together, the ratio A; /A, of the 
perpendicular to the parallel mean free paths remains in the 
range 0.001-0.2. 

In previous studies, a constant solar wind with a typical speed 
of 400km s7}, independent of heliolatitude, and the correspond- 
ing heliospheric magnetic field were usually used in the the- 
oretical analysis of diffusion coefficients and the modulation 
modeling of cosmic rays in the heliosphere. However, a large 
number of spacecraft observations have shown that the real- 
istic solar wind speed varies with the heliolatitude, i.e., with 
low velocity near the ecliptic plane, and with high velocity 
in the polar coronal holes. Therefore, in theoretical investiga- 
tions of heliospheric magnetic field and transport coefficients, 
we should take into account the latitudinal dependence of the 
solar wind speed and its implications for the transport of cos- 
mic rays in the heliosphere, especially in the global modulation 
modeling. 

In this work, we focus on the latitudinally dependent solar 
wind speed and its effects on the heliospheric magnetic field 
and the transport coefficients (parallel, perpendicular, and drift 
scale) of charged energetic particles with adiabatic focusing. 
Specifically, in Section 2, we first present the Ulysses obser- 
vations of the latitudinally dependent solar wind speed from 
1994 September through 1995 July (9/1994-7/1995), known 
as the solar-minimum fast latitude scan. We also provide the- 
oretical models in hyperbolic and piecewise polynomial forms 
to mathematize the latitudinally dependent solar wind speed. 
In Section 3, we investigate the effects and implications of the 
latitudinally dependent solar wind speed on the heliospheric 
magnetic field and the transport coefficients of cosmic rays in 
the presence of adiabatic focusing. A summary of our results 
will be provided in Section 4. 


2. LATITUDINALLY DEPENDENT SOLAR WIND SPEED: 
ULYSSES MEASUREMENTS AND 
THEORETICAL MODELS 


The successful Ulysses mission achieved a number of excit- 
ing new findings regarding the high latitude heliosphere during 
its legendary life. One of the important findings is the latitudi- 
nal dependence of the solar wind speed. Specifically, the solar 
wind is observed to display high speeds in the high heliolati- 
tude coronal holes, low speeds in the coronal streamer belt, and 
intermediate speeds in the transition regions between the coro- 
nal holes and the streamer belt. This complexity was previously 
simplified with a constant and stable solar wind speed through- 
out the heliosphere in the studies of SEP transport and GCR 
modulation. However, on account of the significant difference 
among the solar wind speeds from different coronal sources, 
in the realistic investigations we should consider the latitudinal 
dependence of the solar wind speed. 

We observationally investigate the latitudinally dependent 
solar wind speed measured by the Ulysses spacecraft during 
the solar-minimum fast latitude scan in the time period 9/1994— 
7/1995. In the high heliolatitude regions (both north and south), 
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the solar wind speed is high with a typical value of 800 km s7!; 
however, in the coronal streamer belt near the ecliptic plane, the 
solar wind speed is relatively low with average value of about 
400kms~!. In addition, in the transition regions between the 
coronal holes and the streamer belt, the solar wind speed remains 
in the range 400-800 kms~!. Therefore, the variability of the 
solar wind speed along the heliolatitude is obvious in the global 
view of heliosphere. Apparently, a constant value independent of 
heliolatitude, e.g., the commonly used 400 km s~', is insufficient 
to describe the realistic characteristics of the solar wind speed 
throughout the heliosphere. It is necessary to find a more 
appropriate way to model the latitudinally dependent solar 
wind speed. 

Recently, a hyperbolic model was provided to depict the 
latitudinal dependence of the solar wind speed (e.g., Burger 
& Sello 2005; Heber & Potgieter 2006; Miyake & Yanagita 
2008). This model can be expressed as 


V(0) = Vo (15 + 0.5 tanh [8.0 (6 = A 4 #)]) o (2) 


with Vo = 400kms™! and ¢ = a + B, where a is the tilt 
angle between the rotation and magnetic axes of the Sun, and 
in association with œ is used to determine at which polar 
angle or colatitude (0) the solar wind speed starts to increase 
from 400 km s7! toward 800 km s™!. We choose œ = 2/12 and 
B = x/12 in this work. In Equation (2), the top and bottom 
signs denote the cases of the northern and southern hemispheres 
of the Sun, respectively. 

Alternatively, we provide a piecewise polynomial model to 
quantify the latitudinally dependent solar wind speed. Accord- 
ing to the fact that different solar wind speed categories origi- 
nate from different heliolatitude regions (coronal streamer belt, 
polar coronal hole, and transition region between them) with 
different physical properties, we approximately divide the polar 
angle or colatitude (0) into seven bins, i.e., [0, 40°), [40°, 60°), 
[60°, 80°), [80°, 100°], (100°, 120°], (120°, 140°], (140°, 180°]. 
Within each bin of colatitude 0, we construct a polynomial 
function to model the variation of the solar wind speed with 
the colatitude. The piecewise polynomial function can be 
described as 


800 0 € [0, 40°) 
— 102 + 400 0 € [40°, 60°) 
16? — 800 +3600 0 € [60°, 80°) 
V(6) = 4 400 0 €[80°, 100°] (3) 


56? — 1000 +5400 8 € (100°, 120°] 
—$6? +1400 — 9000 0 € (120°, 140°] 
800 0 € (140°, 180°]. 


In Figure 1, we compare the fitting results of the two theo- 
retical models with the Ulysses observations of the latitudinally 
dependent solar wind speed. The black line indicates the lati- 
tudinally dependent solar wind speed measured by the Ulysses 
spacecraft during the solar-minimum fast latitude scan in the 
time period 9/1994—7/1995. The red and the blue lines indicate 
the hyperbolic and the piecewise polynomial models, respec- 
tively. As one can see, both the hyperbolic and the piecewise 
polynomial models fit to the observations very well. Generally, 
the overall properties of the solar wind speeds originating from 
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Figure 1. Black line denotes the latitudinally dependent solar wind speed measured by the Ulysses spacecraft during the solar-minimum fast latitude scan in the time 
period 9/1994-7/1995. The variability of the solar wind speed along the heliolatitude is obvious in the global view of the heliosphere. In the high heliolatitude regions 
(both north and south), the solar wind speed is high with a typical value of 800kms~!; in the coronal streamer belt near the ecliptic plane, the solar wind speed is 
relatively low with an average value of about 400 km s7!; in the transition regions between the coronal holes and the streamer belt, the solar wind speed remains in 
the range 400-800 km s~!. The red and the blue lines indicate the fitting results of the theoretical models in hyperbolic and piecewise polynomial forms, respectively, 
to the Ulysses observations of the latitudinally dependent solar wind speed. Both the hyperbolic and the piecewise polynomial models fit to the Ulysses observations 


very well. 


different coronal sources have been successfully described. In 
other words, the solar wind is properly characterized with high 
speed in the coronal holes, low speed in the streamer belt, and 
intermediate speed in the transition regions between the coronal 
holes and the streamer belt. 


3. EFFECTS OF LATITUDINALLY DEPENDENT 
SOLAR WIND SPEED 


The classical spiral model of the heliospheric magnetic field 
first suggested by Parker (1958) can be described with 


Bo f. Or sin , 
B=- ê, 7 êp | (4) 


where Bo is a constant chosen to be 3.54 nT AU? in this work 
characterizing the solar magnetic field, V is the solar wind speed 
with a constant value, Q is the angular rotation velocity of the 
Sun, ris the heliocentric radial distance, and 0 is the polar angle. 
The magnetic field magnitude of the standard Parker model can 
be calculated as 


Bo |V? +Q?r? sin? 0 
BC, 6, V)= -3y a (5) 


As we know, the global magnetic field in the heliosphere is not 
constant, and this spatially varying guide magnetic field leads 
to the so-called adiabatic focusing effect. We can derive the 
analytical form of the adiabatic focusing length in the standard 
Parker model as (He & Wan 2012a) 


r(V? + Qr? sin? 0} 


L(r,0, V) = - ; 
r ) V(2V2+2r? sin? 0) 


(6) 


Taking into account the latitudinally dependent solar wind 
speed, Schwadron (2002) and Schwadron & McComas (2003) 


suggested a three-dimensional heliospheric magnetic field 
model, written as 


E Bo (1 rag dV(O@)\., roo » 
a vo do)" Ve? 
(Q — w¢)r sinb , 
V(6) | o 


where wg denotes the rate of magnetic footpoint motion between 
regions of fast and slow solar wind at the Sun, wg denotes 
the differential rotation rate of magnetic field line on the Sun, 
and V(@) denotes the latitudinally dependent solar wind speed, 
which can be empirically expressed with Equation (2) or (3). 
The magnitude of the Schwadron field can be expressed as 


2 2 . 2 
rag dV(0) ro (Q—a¢ )r sin@ 
(1 ZO do ) + (32) +( VO) i 
(8) 


By assuming that wọ ~ Q/4 and rQ/V(@) ~ rp, and setting 
@¢ = 0 in Equation (7), Burger & Sello (2005) obtained a 
simplified two-dimensional model of the Schwadron field as 


_ Bo rp dV(@)\,  Qrsind, 
B= an [0 4V (0) do Ja |: i 


vo 
where r, is a dimensionless constant chosen to be different 
values characterizing different dependence. In this simplified 
two-dimensional version, the Bọ component is also ignored. 
For simplification and convenience, as a first result of this paper 
we adopt the simplified two-dimensional magnetic field model 
used by Burger & Sello (2005) and Miyake & Yanagita (2008). 
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Figure 2. Effect of the latitudinally dependent solar wind speed shown in Figure | on the heliospheric magnetic field magnitude at 1 AU. The black dotted curve 
indicates the magnetic field magnitude of the standard Parker model, while the black (blue) dash-dotted, dashed, and solid curves indicate the magnetic magnitudes 
of the modified two-dimensional (three-dimensional) field model with the dimensionless constant rp = 1, rp = 2, and rp = 4, respectively. 


We will consider later a three-dimensional heliospheric mag- 
netic field taking into account the latitudinally dependent solar 
wind speed. Accordingly, the latitudinally dependent magni- 
tude of the modified two-dimensional heliospheric magnetic 
field could be obtained as 


_ Bo rp dV(@)\> (Qrsind\? 
A= z (i 4V0) dé ) (Se). go) 


We can also derive the latitudinally dependent focusing length 
of the modified two-dimensional heliospheric magnetic field 


4 dé 4V8) dé 


s 2 2 
| (vo-44) 1.0% si 0| v2@)90%7 sin? /(1- Fay Z2) 
Lr, 0) = 


2 
vofa vo-4 KOY 4 0%? sin? e] 


a1) 
Recently, Hitge & Burger (2010) used a three-dimensional he- 
liospheric magnetic field to investigate the cosmic ray modula- 
tion with a latitude-dependent solar wind speed. For comparison 
with the simplified two-dimensional model (Equation (9)), we 
similarly consider a simplified three-dimensional heliospheric 
field model, which is based on the fully three-dimensional 
Schwadron field, as 


a paii 


BE Bo i Tp 8 Or sind , 
~ 72 4Vv(0) do 4% 


ZON 
(12) 


We note that in the ecliptic, spacecraft observations have shown a 
negligible meridional component of heliospheric magnetic field. 
In order to conform to the observations and other models near the 
ecliptic, we recast the simplified three-dimensional heliospheric 
magnetic field model (Equation (12)) into the normalized 


form as 


Bo rp dV(0)\n lpn Qr sinb A 
72 [( W@ da Je, 7&9 V(0) ês] 


at high latitudes, 
B= (13) 


Bo rp dV(@\ A^ Or sin A 
re (i W® dé )é VO) ês | 


near the ecliptic. 


For the simplified three-dimensional heliospheric magnetic 
field at high heliolatitudes, we can also obtain its latitudinally 
dependent magnitude 


EE (1 lp ON" (By + (Seat), 
r2 4V(0) do 4 V0) 
(14) 


and its latitudinally dependent focusing length 


L(r, 0) 


2 2 EEE P a 
ri({1 rp_ dV(0) +242 rf sint 0 
WO do 16 VZO) 
2(1——"P_ dV) 2 TP , 92r? sin? 6 1 _"P_ av qpr? sin 26 i 
4V(0) dé 8 v2) 4V(0) do 8V2(0) 


(15) 


We note that the approaches and the main conclusions presented 
here can be easily extended to the fully three-dimensional 
version of the Schwadron field model (Equation (7)) and other 
heliospheric magnetic field models. 

In the following, we will directly investigate the effects of 
the latitudinally dependent solar wind speed on the magnetic 
magnitude and the adiabatic focusing length by comparing 
the two cases of the standard Parker model and the modi- 
fied two-dimensional or three-dimensional heliospheric mag- 
netic field mentioned above. Figure 2 presents the magnetic 
field magnitude at 1 AU as a function of the colatitude 0 
ranging from 0° to 180°. The black dotted line indicates the 
magnetic field magnitude of the standard Parker model, while 
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Figure 3. Effect of the latitudinally dependent solar wind speed shown in Figure | on the adiabatic focusing length of the heliospheric magnetic field at 1 AU. The 
black dotted curve indicates the adiabatic focusing length of the standard Parker model, while the black (blue) dash-dotted, dashed, and solid curves indicate the 
adiabatic focusing lengths of the modified two-dimensional (three-dimensional) magnetic field model with the dimensionless constant rp = 1, rp = 2, and rp = 4, 


respectively. 


the black dash-dotted, dashed, and solid curves indicate the 
magnetic magnitudes of the modified two-dimensional field 
model with the dimensionless constant rp, = 1, rp = 2, and 
rp = 4, respectively, and the blue dash-dotted, dashed, and 
solid curves indicate the magnetic magnitudes of the modified 
three-dimensional field model with the dimensionless constant 
rp = l, rp = 2, andr, = 4, respectively. In the figure one can 
see that there exists a significant difference between the standard 
Parker field and the modified magnetic fields with latitudinally 
dependent solar wind speed, especially within the two transi- 
tion regions between the coronal holes and the streamer belt. 
This indicates that the effect of the latitudinally dependent 
solar wind speed on the heliospheric magnetic field magnitude 
is very considerable. 

Figure 3 presents the adiabatic focusing length of the helio- 
spheric magnetic field at 1 AU as a function of the colatitude 
0 ranging from 0° to 180°. The black dotted curve indicates 
the adiabatic focusing length of the standard Parker model, 
while the black dash-dotted, dashed, and solid curves indicate 
the adiabatic focusing lengths of the modified two-dimensional 
magnetic field model with the dimensionless constant r, = 1, 
rp = 2, andr, = 4, respectively, and the blue dash-dotted, 
dashed, and solid curves indicate the adiabatic focusing lengths 
of the modified three-dimensional magnetic field model with 
the dimensionless constant r, = 1, rp = 2, andr, = 4, re- 
spectively. We can see that the latitudinally dependent solar 
wind speed significantly affects the adiabatic focusing length 
of the heliospheric magnetic field. We also notice that in the 
modified two-dimensional magnetic field model, the choice of 
the value of r, is not very important. However, in the modified 
three-dimensional model, the choice of the value of r, is still 
important when we consider the adiabatic focusing length in the 
high heliolatitude regions. 

The diffusion tensor, consisting of parallel diffusion, perpen- 
dicular diffusion, and drift coefficients, is an essential parameter 
to study the transport and modulation of charged energetic par- 
ticles in the turbulent interplanetary and interstellar magnetic 
field. In what follows, we will investigate the effects of the 


latitudinally dependent solar wind speed on the diffusion co- 
efficients of cosmic rays in the presence of adiabatic focusing. 
First, we employ the description of the parallel mean free path 
of cosmic rays with adiabatic focusing (see He & Schlickeiser 
2014, and references therein): 


oy L 
Ay =3L [cot ($) — | , (16) 
L Ào 


where Ao is the parallel mean free path in the absence of adiabatic 
focusing, Lis the adiabatic focusing length of the magnetic field. 
By combining Equation (16) with Equation (15), we obtain 
the parallel mean free path of cosmic rays in the presence of 
adiabatic focusing with the effect of the latitudinally dependent 
solar wind speed. We can also investigate the effects of the 
latitudinally dependent solar wind speed on the perpendicular 
mean free path and the drift scale of cosmic rays in the 
heliosphere. We use an analytical form of the perpendicular 
mean free path as (e.g., Burger et al. 2008; He & Wan 2012b) 


2 
—~1 _T(44+1). 3B2,]3 
AL = at Wha ( Di 2D GADS, (17) 
2q + ) 


where a is a constant chosen to be 1/ J3, bp and ô Be are the 
bendover scale and the variance of the two-dimensional compo- 
nent of the turbulence, respectively. In this work, we employ the 
slab/two-dimensional composite magnetic turbulence model, in 
which the turbulence energy in the two-dimensional and slab 
modes is in the ratio of 80:20. We also utilize the typical so- 
lar wind conditions at 1 AU: Vo = 400km s7}, the turbulence 
strength 6B/B = 0.4, the slab bendover scale lsa = 0.03 AU, 
the two-dimensional bendover scale lap = 0.1/5 = 0.003 AU, 
and the inertial range spectral index q = 5/3. The so-called 
classical scattering result of the drift scale of cosmic rays is 
written as (e.g., Burger & Visser 2010; Tautz & Shalchi 2012) 


ATRL 
= 32 2° 
Ài + Ry 


da (18) 
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Figure 4. Effects of the latitudinally dependent solar wind speed shown in Figure 1 on the parallel mean free path (top panel), perpendicular mean free path (middle 
panel), and drift scale (bottom panel) of cosmic ray protons with energy of 100 MeV at 1 AU. In each panel, the black dotted curve indicates the corresponding 
diffusion coefficient associated with the standard Parker model, while the black (blue) dash-dotted, dashed, and solid curves indicate the corresponding diffusion 
coefficients associated with the modified two-dimensional (three-dimensional) magnetic field model with the dimensionless constant rp = 1, rp = 2, and rp = 4, 


respectively. 


where Rz is the unperturbed Larmor radius of the cosmic ray 
particle. 

Figure 4 presents the parallel mean free path, perpendicu- 
lar mean free path, and drift scale of cosmic ray protons with 
energy of 100 MeV at 1 AU, as a function of the colatitude 
0, from top to bottom panels, respectively. In each panel, the 
black dotted curve indicates the corresponding diffusion co- 
efficient associated with the standard Parker model, while the 
black dash-dotted, dashed, and solid curves indicate the cor- 
responding diffusion coefficients associated with the modified 
two-dimensional magnetic field model with the dimensionless 
constant r, = 1, rp = 2, andr, = 4, respectively, and the 
blue dash-dotted, dashed, and solid curves indicate the cor- 
responding diffusion coefficients associated with the modified 
three-dimensional magnetic field model with the dimensionless 
constant r, = l, rp = 2, andr, = 4, respectively. From the 
figure one can see that the diffusion coefficients derived from 
the standard Parker field gradually change with the colatitude 
0 from 0° to 180°. However, the diffusion coefficients derived 
from the modified heliospheric magnetic fields with the effects 
of the latitudinally dependent solar wind speed vary more sig- 
nificantly through the colatitude range of 0°-180°, especially 
within the two transition regions between the coronal holes and 
the streamer belt, where the diffusion coefficients of the modi- 
fied magnetic fields are much smaller than those of the standard 
Parker field. This implies that the latitudinally dependent so- 
lar wind speed plays a very important role in the diffusion and 
transport processes of cosmic rays in the global heliosphere. 
Therefore, when we investigate the transport and modulation of 
cosmic rays in the global view of the heliosphere, the effects 
of the latitudinally dependent solar wind speed should be taken 
into account. 


4. SUMMARY AND CONCLUSION 


The Ulysses spacecraft has shown that the solar wind speed 
varies with the heliolatitude, i.e., high speed in the high heliolat- 
itude coronal holes, low speed in the coronal streamer belt, and 


intermediate speed in the transition regions between the coronal 
holes and the streamer belt. Therefore, a constant solar wind 
speed independent of heliolatitude is oversimple for describing 
the realistic conditions. It is necessary to theoretically model 
the latitudinally dependent solar wind speed and systematically 
investigate its effects on the heliospheric magnetic field and the 
transport coefficients of cosmic rays in the presence of adiabatic 
focusing. In this work, we present theoretical models in hyper- 
bolic and piecewise polynomial forms to model the latitudinally 
dependent solar wind speed measured by the Ulysses spacecraft 
in the time period 9/1994—7/1995. Subsequently, we further 
investigate the effects of the latitudinally dependent solar wind 
speed on the heliospheric magnetic field and the diffusion coef- 
ficients (parallel, perpendicular, and drift scale) of cosmic rays 
with adiabatic focusing. We find that the latitudinally dependent 
solar wind speed has significant influences on the magnitude of 
the heliospheric magnetic field, especially within the two tran- 
sition regions between the coronal holes and the streamer belt. 
In addition, the latitudinally dependent solar wind speed plays 
a very important role in the diffusion and transport processes 
of cosmic rays in the presence of adiabatic focusing in the he- 
liosphere, especially within the two transition regions between 
the coronal holes and the streamer belt. Therefore, we suggest 
that in the three-dimensional or global heliospheric modeling 
studies of cosmic rays, the effects of the latitudinally dependent 
solar wind speed should be taken into account. In the theoret- 
ical analysis of cosmic ray diffusion and transport, we should 
also employ the modified descriptions of the heliospheric mag- 
netic field and the diffusion coefficients of cosmic rays with the 
effects of the latitudinally dependent solar wind speed. 

In most large SEP events, the SEP flux increases are usually 
observed concurrently by in-ecliptic near-Earth spacecraft and 
the Ulysses spacecraft at high heliolatitudes (e.g., McKibben 
et al. 2003; Lario et al. 2003; Dalla et al. 2003; Sanderson 
et al. 2003). Some simulation works have been presented to 
explain this phenomenon (e.g., He et al. 2011). In the studies of 
SEP transport and GCR modulation in the future, we will use 
the modified expressions of the diffusion coefficients (parallel, 
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perpendicular, and drift scale) of charged energetic particles 
in the presence of adiabatic focusing and with the effects of 
the latitudinally dependent solar wind speed, provided in the 
present paper. In addition, we note that the theta component of 
the modified heliospheric magnetic field structure is azimuthally 
phase dependent, which means that it has the potential to create 
significant asymmetries. In the future work, we will also take 
into consideration these complications and their effects on 
the drifts and parallel propagation of charged particles in the 
heliosphere. 
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